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ABSTRACT 

We present phase-resolved spectroscopy and photometry of the optical counterpart to 
the X-ray binary IGR J19308-I-0530. Ellipsoidal modulations in the lightcurve show 
that the F-type companion star in the system is Roche-lobe filling. The optical spectra 
are dominated by absorption features from the donor star, with ~ 10-20% disc contri- 
bution to the optical continuum. We measure an orbital period of 14.662±0.001 hours, 
a radial velocity semi-amplitude for the companion star of K2 =91.4±1.4kms~^and a 
rotational broadening of v sin i = 108.9±0.6 km s~^ . From K2 and v sin i, given that the 
donor star is filling its Roche lobe, we derive a mass ratio of q =M2/Mi = 1.78±0.04, 
which is typically considered to be too large for stable Roche-lobe overflow. Our ob- 
servations support an inclination of ^50 degrees. The accretor in IGR J19308-1-0530 
is most likely a white dwarf, although a neutron star cannot entirely be excluded. 

Key v^rords: stars: individual (IGR J19308-1-0530) — accretion: accretion discs — 
stars: binaries — X-rays: binaries 



1 INTRODUCTION 

Intermediate-mass X-ray binaries (IMXBs) are binary sys- 
tems where a compact object - black hole, neutron star (NS) 
or white dwarf (WD) - is accreting matter from a compan- 
ion star of spectral type A or F. IMXBs are rarely observed 
(see, e.g., the catalogue from Liu et a l.'l2007" and the 2012 
version of the RK catalogue [Ritter fc Kolb.200jy] ). The ma- 
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jority of accreting WDs in binaries, in fact, belongs to the 
class of cataclysmic variables (CVs), which have late type 
secondaries (M2 < Mi) and -P < 6 hours (Knigge et al. 
|2011[ ). NSs or BHs, instead, are typically observed in X- 
ray binaries (XRBs) hosting either a massive 0-B donor 
star (high-mass XRBs, HMXBs) driving accretion via stel- 
lar wind, or a late M or K dwarf secondary star (low-mass 
XRBs, LMXBs) accreting via Roche lobe overflow. The rea- 
son for the observed rarity of IMXBs especially among NS 
and WD systems is that, when the companion is more mas- 
sive than the accretor but not massive enough to have strong 
winds, wind accretion proceeds at a very low rate and Roche- 
lobe accretion is thought to be unstable. For NSs and WDs 
in IMXBs mass flows from the more massive to the lighter 
star and angular momentum conservation shrinks the orbit, 
leading to enhanced mass transfer. The bright X-ray binary 
phase is therefore intense and short-lived, causing an obser- 



vational bias towards LMXBs, CVs and HMXBs (Tauris & 
[van den Heuvel|2006[ ). Nevertheless, IMXBs could be a large 
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fraction of the XRB population and have an important role 



in understanding their evolution ( Podsiadlowski et al. 2001 1 



Cyg X-2 and Her X-1 are though to have started as IMXBs, 
even though the measured mass ratio is currently < 1. 



IGR J 1930 8+0530 was discovered by INTEGRAL ([Bird 
et al.||2006| ) and observed by Swift ( [Rodriguez et al.||2008| ). 

An association with the star TYC 486-295-1, classified eis 
an F8 star in the survey by 'McCuskey (19491, was made 



using the Swift position ( ^Rodriguez et al. 20081. This was 
confirmed using an accurate Chandra position of the X- 



ray source (Ratti et al. 20101. Considering typical param- 



eters of an F8 star, Rodriguez et al. ( 2008 1 suggested 



IGR J19308-f0530 to be a L/IMXB in quiescence or a CV 
at a distance of < 1 kpc. 

Here, we present phase-resolved optical spectroscopy and 
photometry of IGR J19308-I-0530, in order to measure the 
orbital period P, the radial velocity semi-amplitude K2 and 
the projected rotational velocity v sin i of the companion 
star, and the system inclination i. In a Roche lobe filling 
system, K2 and v sin i allow us to infer the ratio q —M2 /Mi 
between the mass of the secondary and the primary star in 
the system (Wade & Hornc 1988, see also jGray^ 1992J and, 
knowing P and i, to solve the system mass function. 



2 OBSERVATIONS AND DATA REDUCTION 

In total twenty-two high-resolution spectra of 
IGR J19308-I-0530 were collected. Observations were 
made on two nights in 2010 Mar., one night in 2010 
Apr. and 9 nights in 2010 Jun. using the fiber spectrograph 
High Efficiency and Resolution Mercator Echelle Spec- 
trograph (HERMES), mounted at the Mercator telescope 
in La Palma ( [Raskin et al. 20111. The typical exposure 
time was 1200 seconds. The spectra have a dispersion of 
0.027 A/pixel at 5000 A and cover the wavelength range 
3770-7230 A. The fiber aperture is 2'.'5 on the sky, but 
the presence of a slicer mimics a narrow slit providing 
a resolution of ~85000 irrespective of the seeing. The 
template star HD185395, of spectral type F4 V (later found 
to be the closest match to the spectral type of the target, 
see below) was observed with the same settings for 360 
seconds on 2009 Aug. 5. The extraction of the spectra was 
performed through the dedicated automated data reduction 
pipeline HermesDRS. For each spectrum we selected two 
regions for the analysis, one covering the H7 and H/3 lines 
(4280-5250 A) and one around the Ha line (5950-6700 A), 
which we will refer to as SI and S2, respectively. These 
regions were selected as they are rich in stellar lines, with 
little contamination from interstellar features. A good 
fit of the contirmum was achieved in each region with a 
polynomial function of order 9. We normalized the spectra 
dividing by the polynomial fit. 

We also performed time-resolved photometry of 
IGR J19308-I-0530, with the 80-cm IAC80 telescope at 
the Observatorio del Teide in Tenerife equipped with the 
CAMELOT CCD imager. The observations were obtained 
during part of three nights between 2010 Jul. 30 and Aug. 02, 
by cycling through the Sloan g' , r' , i! and ^'-band filters. 
Three consecutive exposures were obtained in each filter, 
with integration times ranging from 6 to 24 sec depend- 
ing on the filter and seeing conditions. Standard stars were 
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Figure 1. Optimal subtraction of the spectrum of the F4V stellar 
template HD185395 from the spectrum of IGR J19308-(-0530. The 
image shows the part of spectral region S2 used to measure v sin i. 
The interstellar feature at '~6280 A was masked. The template 
and target spectrum are offset by 0.4 and 0.2 in y. 



not taken due to non-photometric weather. After debias- 
ing and flat-fielding the images using dome fiat-field obser- 
vations (with standard routines in IRAf) the instrumental 
magnitudes of IGR J 19308-1-0530 and four comparison stars 
were computed by means of aperture photometry. Differen- 
tial lightcurves were then obtained for IGR J19308-I-0530 
with respect to the comparison star TYC 486-968-1. We 
also extracted lightcurves for TYC 486-968-1 using the other 
comparison stars. No signiflcant variability was detected 
in the i'-band (r.m.s.< 0.006 mag), whereas for the other 
bands a larger scatter was observed due to weather condi- 
tions (clouds and Calima). The maximum departure from 
the mean value was of 0.1 mag. Given that the ellipsoidal 
modulation in IGR J19308-I-0530 is small in amplitude we 
decided to model only the i'-band lightcurves. 



3 ANALYSIS AND RESULTS 
3.1 Spectroscopy 

The spectra are dominated by absorption features from the 
secondary star (??). Although no emission line is directly 
visible, residual emission in the Balmer lines appears when 
subtracting the spectra one from another, after correcting 
for the orbital shift of the lines. The emission component 
is variable in intensity and wavelength and slightly shifted 
with respect to the absorption line. 

Using the package molly developed by T. Marsh, we 
measured the orbital velocity of the companion star in 
IGR J19308-f0530 by cross-correlating the spectra of the 
target with that of the template star. As the absorption lines 
are rotationally broadened in the target spectrum, a broad- 
ening of 100 kms~^was applied to the template (see below 
on rotational broadening), improving the cross-correlation. 
The Balmer lines and interstellar features were masked. 
We performed a fit of the velocities versus time with a 
sine function, with K2, the systemic radial velocity 7, To, 
and P as free parameters. To was constrained to be near 
the middle of the time span over which the observations 
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Figure 2. Radial velocity curve from the absorption lines in 
the IGR J19308+0530 spectra (region SI) and its best fitting 
sinusoid (solid line). For comparison we plot circular orbits with 
K2 =112 km s^i (dotted Une) and K2 =140 km s"^ (dashed line), 
providing g = 1.4 and q = 1 respectively for the observed vsini. 



were taken, and such that phase is at the inferior con- 
junction of the companion star. The best-fitting sinusoid 
provided a of 57.6 for the region SI, and 41.76 for S2 
(19 d.o.f.). In both cases, the uncertainties on the param- 
eters were estimated assuming that the sinusoidal model 
was correct, and we scaled the errors on the velocities to 
reach a reduced of ~1. The values of P, Tq and 7 mea- 
sured in SI and S2 are consistent at the Icr level, K2 is 
consistent at the 2 a level. The error-weighted average of 
the parameters gives P = 0.61092 ± 0.00003 days. To = 
2455330.8169 ± 0.0023 HJD/UTC, K2 = 91.4 ± 1.4 kms"^ 
and 7 = —18.5 ± 0.9 kms~^. Figure [2] shows the radial ve- 
locity curve (rvc) folded on the above period. The value 
of 7 is in the reference frame of the template star used 
for the cross-corr elation, whose systemic radial velocity is 
-28 ± 0.9 kms~^( |wiison||l953[ ). The systemic radial veloc- 
ity of IGR J19308-h0530 is therefore -46.5 ± 1.2 kms'V 

We obtained a set of high signal-to-noise template 
spectra from the UVES Paranal Observatory Project 
(UVESPOP, Bag nulo et al.f 2003) of A, F and early G stars 
of luminosity class V, III and IV. We subtracted each tem- 
plate spectrum from the Doppler-corrected average of the 
IGR J19308-I-0530 spectra between orbital phase 0.9 and 
0.1. The reason for choosing this range in phase is that the 
oblate shape of the Roche-lobe filling companion star and 
the possible presence of irradiation from the compact ob- 
ject could cause asymmetries in the line profiles, which are 
minimised close to phase 0. We performed a test on the 
residuals of the subtraction: the template resulting in the 
minimum provides our best estimate for the source spec- 
tral type. In particular, we adopted the optimal subtraction 
procedure implemented in molly, where the templates are 
multiplied by a factor 0</< 1 before the subtraction, rep- 
resenting the fractional contribution of light from the sec- 
ondary star (1 minus the disc veiling). The factor / is found 
by minimizing the difference between the residuals and a 
smoothed version of itself. Before doing the subtraction, the 
UVESPOP spectra were shifted to the rest frame of the 
average target spectrum, and degraded to match the sam- 
pling and line broadening of the latter. 

The procedure favors an F4-F6 V companion star, with a 
disc veiling of 10 — 20%. The same spectral type and veil- 
ing are obtained when considering an average of the target 



spectra around phase 0.5, suggesting little irradiation on the 
inner face of the companion star. 

To measure vsini, we compared the spectrum of the 
template star with the Doppler-corrected average of the 
IGR J19308-I-0530 spectra between phase 0.9 and 0.1. 
The observed full-width at half-maximum (FWHM) of the 
absorption lines in the target spectra is determined by 
the intrinsic line width (expected to be dominated by 

V sin i) , broadened by the instrumental resolution profile and 
smeared by the motion of the companion star during the in- 
tegration time of one observation. In order to account for 
the smearing, we made as many copies of the template spec- 
trum as the number of target spectra we used for the av- 
erage and we artificially smeared each copy of the template 
by 2-nTK cos(2Ti(f)) / P , where T is the duration of one ex- 
posure on IGR J19308-I-0530 and the phase of one of the 
IGR J19308-I-0530 spectra we averaged. After that, we av- 
eraged the smeared template and broadened the resulting 
spectrum with different values of vsini. For each wsini, we 
performed an optimal subtraction of the broadened template 
from the averaged spectrum of IGR J19308-I-0530: again the 
broadening which gives the minimum provides a measure 
of the actual vsvcii (Figure [T]). A value of 0.5 was assumed 
for the limb darkening. We masked interstellar features and 
the lines from the Balmer series. 

In order to estimate the uncertainty on «sinj, we in- 
cluded this procedure in a Monte Carlo simulation, following 
Steeghs fc Jonker] ( |2007[ ). We copied each target spectrum 
500 times, using a bootstrap technique where the input spec- 
trum is resampled by randomly selecting data points from 
it. The bootstrapping maintains the total number of data 
points in the spectrum. For each bootstrap copy, one value 
of «sinj is measured as described above. The distribution 
of V sin i obtained from the 500 copies is well described by a 
Gaussian, whose mean and r.m.s. provides the best-fit i; sin i 
and its Icr error. As template and target spectra are acquired 
with the same instrument, the instrumental resolution pro- 
file is not affecting our measurement. 

The weighted average of the results from SI and S2, con- 
sistent at the Icr level, is i;sini = 108.9 ± 0.6kms~^. With 

V sin i and K2 , we calculated the system mass ratio q — 



M2/M1 from the relation 



_ ^ = (l + ^)a 

( Home et al.|1986| ), obtaining q =1.78±0.04. 



6ij2/3+ln{l + ql/3) 



3.2 Ellipsoidal modulation and system inclination 

Figure [3] shows the i'-band hghtcurve for IGR J19308-I-0530 
obtained by phase folding the CAMELOT data on the 
ephemeris determined in Section |3.1[ The lightcurve dis- 
plays the typical signature of ellipsoidal variation, with 
two unequal minima, but in addition asymmetric maxima 
( 0'Connell|l95T1 and [Wilsey fc Beaky|2009 |. 
We modeled the i'-band lightcurve using the XRbinary pro- 
gram written by E.L. Robinson. A reasonable fit to the data 
(reduced ~ 7, 403 d.o.f.) is obtained with a model assum- 
ing an F4V secondary star and including 30% disc contri- 
bution to the total light plus a disc hot-spot at phase 0.75. 
The disc veiling in the model is larger than observed in the 
spectra, but variability is possible as the photometry was 
performed one month after the last spectrum was acquired. 
Similar models with different assumptions about the disc 
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Figure 3. j'-band CAMELOT lightcurvc, folded using the 
ephemeris from the radial velocity data. The colours indicate the 
three observing nights (Jul. 30 in red, Jul. 31 in green and Aug. 1 
in blue). The drawn line shows a fit with ellipsoidal modulations 
plus a disc with a bright disc spot. 



properties (disc radius, height and temperature profile) also 
give reasonable fits. As we have no indications to single out 
one preferred set of parameters, we do not provide a formal 
uncertainty on i, but an indicative value of ~52 °as a guide. 



3.3 Upper limits in the radio waveband 

We searched for a radio counterpart to IGR J19308+0530 
with a Karl G. Jansky Very Large Array observation at 4.6 
and 7.9 GHz (120 MHz bandwidth) (proposal ID lOB-238), 
taken on 2010 Aug. 19 with the array in its most compact 
D-configuration. The on-source time was 33 minutes and 
the data were reduced according to the standard procedures 



within the Common Astronomy Software Application ( Mc- 



MuUin et al. 20071 software package, using the calibrator 



3C 286 to set the amplitude scale and J1922+1530 to cali- 
brate the amplitude and phase gains for the target source. 
IGR J19308+0530 was not detected at either frequency. A 
3-sigma upper limit to the source flux of 54 /xjy/beam was 
derived 

To search for radio pulsations from a potential pulsar in 
this system, we used the Westerbork Synthesis Radio Tele- 
scope and the PuMa2 backend ( Karuppusamy et al.||20d8 |. 
We observed for 1 hr from both 310-380 MHz (Obs. ID 
11205182, 29 Aug. 2012) and from 1300-1460 MHz (Obs. 
ID 11205210, 31 Aug. 2012). The predicted dispersion mea- 
sure (DM) for a 300 pc distance along this line-of-sight (see 
Discussion) is only ~ 4pc cm~^ ( Cordes fc Lazio|2002 l. For 
each data set, we used the PRESTO software suite jRansom 



2001 1 to search a set of trial DMs up to 30 pc cm . No ob- 



vious radio pulsar signal was detected after an acceleration 
search. From these observations, we can place conservative 
flux density limits of 6*350 < 0.9 mjy and 51400 < 0.4 mjy 
for any pulsar present in the system, assuming it is beamed 
towards us. 



4 DISCUSSION AND CONCLUSION 

We performed a dynamical study of the system 
IGR J 19308+0530 through optical spectroscopy and 
photometry. The optical spectra are dominated by the 
companion star, with no evidence of irradiation, no 
emission features visible from the accretion flow besides 
a partial flUing in of the Balmer lines and 10 — 20% 
disc contribution to the continuum. The secondary star 
is most likely of spectral type F4-6V. Ellipsoidal mod- 
ulations are detected on the ~14.6 hour orbit. From 
phase-resolved spectroscopy we measure an extreme value 
for the binary mass ratio of q =1.78±0.04. The lightcurve 
modeling provides a reasonable fit to the data with a 
disc+hot-spot model at i ~ 52°. Solving the mass function 
/(M2) : Mi^f^ = =0.03 A/0 with this inclination, 

we obtain the following indicative masses: Mi ~0.8Mq and 
M2 ~1.4 Mq. The masses are consistent with a WD accretor 
and an F4V donor (typical mass of ~1.37A/0, Mamajek's 
list 20110. If the inclination is lower, which we cannot 
exclude based on these data, the masses will increase, 
allowing a scenario with a NS primary if i < 45°. However, 
in this case the companion star would be over-massive for 
the spectral type, which is unusual for XRB and CVs. 

Assuming an F4V mass donor with a radius equal 
to that of the Roche lobe for our best estimated masses, 
the magnitude of the companion is Mv ~ 2.7 in the vi- 
sual band (half a magnitude brighter than for a typical 
F4V star, Mamajek 2011). Comparing with the apparent 
magnitude of IGR J19308+0530 {mv = 10.95 ± 0.13, con- 
verted from nivT in the Tycho catalogue), we estimate a 
distance range of 300 — 45 pc, for an extinction b etween 
Nh = (0-2.6) X 10^^ cm"^ (iDickey fc Lockman|l990l where 



Nh is converted into Ay following Giiver fc Ozel|20'09 1 



Combining the systemic radial velocity with the source 
proper motion reported in the UCAC3 catalogue (?), 
we computed the Galactic space velocity components 
of IGR J 1 9 308-1-0 530 using the method of [Johnson fc] 
Soderblom (19871. Assuming that the Local Standard 



of Rest ( LSR) participates in the Galactic rotation at 
254kms~^ |Reid et al.|[2009| ) and a distance of 375 ± 75 pc, 
the derived peculiar velocity is 45.3±2.9kms^^. This rules 
out a large asymmetric kick from a supernova. 

The X-ray lumin osity of IGR J193 08-I-0530 measured 



by Chandra in 2007 ( [Ratti et al.||2010[ ) is 5 x lO^"* - 4 x 
lO^^ergs"^, for distance in the above range of 300 — 450pc 
and the corresponding Nh- In the same way, UV Swift ob- 
servations provide 4 x 10^" - 1 x lO^^ergs"^ A"^ at 2500 A. 
A scenario where IGR J 19308-1-0530 is not in full contact 
and the X-rays are from coronal activity of the compan- 
ion seems unlikely. First, coronal activity usually does not 
produce prolonged high energy emission, while the source 
was discovered from a stack of INTEGRAL observations. 
Second, coronal activity could explain the observed X-ray 
luminosity only for the earliest spectral type allowed by 
the observations combined with very little Nh, as for F 
stars the ratio between the X-rays and bolometric flux is 
log{Lx/Ltoi) < -4.6 for t;sini ~ 100 (|Walter|[T983l). Fi- 



^ http://www.pas.rochester.edu/ ^emamajek/ 
EEM_dwarf_UBVIJHK_colors_Tefr.dat 
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nally, the spectra indicate disc contribution to the contin- 
uum, and the emission detected in the Balmer hues shows a 
slight, variable velocity offset with respect to the radial ve- 
locity curve of the companion star which is not expected in 
case they originate from coronal activity. The emission line 
component also seems weak for an highly active star. A hot 
WD of --^60000 K alone could account for the UV emission, 
but not for Lx unless thermonuclear burning is happening 
on its surface. The few intermediate-mass CVs with long or- 
bital period that are known do appear as super soft sources 
(SSSs), showing soft X-ray spectra possibly due to stable hy- 
drogen burning on the WD (e.g. [Ka habka 2006|. However, 
with a luminosity of 10^'' -10^* ergs"', SSSs are much more 
luminous than what we observe from IGR J19308-I-0530. As 
the wind mass loss expect ed for the secondary spe ctral type 
is low (~ 10""Mq yr' MCranmer fc Saar||201l[ ) wind ac- 
cretion can not account for the observed X-ray luminosity. 
We conclude that IGR J19308+0530 is most likely in con- 
tact, consistent with the ellipsoidal modulations observed in 
the lightcurve. This makes the source particularly interest- 
ing among the scarcely populated class of IMXBs since it 
shows Roche lobe overflow at a low accretion rate m. This 
is unusual, as Roche lobe accretion with an intermediate 
mass companion and such a large mass ratio is typically 
considered unstable, with an intense and short-lived accre- 
tion phase (e.g., Tauris fc van den Heuvel|20 06|). Donor stars 
of mass around 1.4 have thin or non-existent convective 
envelopes ( Verbunt fc Zwaan||1981 1 which implies that the 
instability does not proceed on a dynamical timescale. While 
a convective envelope star would react to mass loss by ex- 
panding, a star with a mostly radiative envelope will shrink, 
slowing down the mass transfer. However, with a mass ra- 
tio of 1.8 the radius of the Roche lobe will reduce faster 
in response to mass transfer than the radius of the com- 



panion, triggering a thermal instability (Tauris fc van den 
|Heuvel|[20 06 ) . Using the binary stellar-evolution code orig- 
inally developed by Eggleton ( Yakut fc Eggleton||2005 and 
references therein), we modeled evolutionary tracks that al- 
low for periods of stable Roche-lobe accretion onto a WD 
with masses and orbital periods consistent with our find- 
ings, although at higher m values than implied by the Lx 
(see also the appendix in [Podsiadlowski et al.|[2003[ ). The 
low Lx can be explained if we are observing a short lived 
phase of low m, or if the mass transfer is non-conservative 
(a fully non-conservative scenario requires a mass loss rate 
from the system of 10"^" — lO~^M0/yr). There is currently 
no direct evidence in favour of or against the presence of 
mass loss from the system. 

From an observational point of view, an overestimate 
of q could be due to uncertainties on the limb darkening. 
However, even assuming a limb darkening of (instead of 
0.5) reduces usini by only a few kms~^, still providing a 
high g of ~ 1.6. 
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